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** Abstract : 

-Calculating the maximum value Pmax that causes the fracture propagation while applying the internal 
static pressure . 

-Evaluating the time-limit for applying the continous / discrete form of internal dynamic pressure for the 
round metal vessel under conditions given by series expansion and curve -fitting method . 

** Subjects: Fracture Mechanics , The stress intensity factor KI . 



Introduction 

This worksheet demonstrates Maple's capabilities in estimating the maximum pressure and the time- 
limit for applying static or dynamic pressure on the surface of a round metal vessel which has an internal 
semi-elliptical surface crack . 

All rights reserved. Copying or transmitting of this material without the permission of the authors is not 
allowed . 



1. Modeling the problem 



We consider the round metal vessel with parameters given : 

a : one-half crack depth (mm) ; B : thickness of vessel cm) ; c : one-half crack length (mm) ; d : internal 
diameter(m) ; P : internal pressure (MPa) . 



Modeling the sizes of the vessel and crack 




The vessel with an internal semi- elliptical crack 
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2. Technical basis and calculation formulas 

Depending on the principle of superposition , the total stress intensity factor will be determined by 



K 



Itotal 



K + K 

I Hoop I Pressure 



After calculating , K can be rewritten in the form : Kj. . , 

5 ' Itotal ^'"*"' 



aP -J an 

pressure ' 



a : correction/actor 



1.12 



/§ 



shape 



diameter 
thickness 



+ 1 



Use the formulas of hoop stress : a 



1 press diameter 



hoop 



B 



and 



thickness 



the shape factor Q^^^^^ = — n 



1 a 



2 \ 
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, we find the value of correction factor 



a , consequently it gives the total stress intensity K . Because the critical crack length is 

Itotal 



defined by : a 



crack 






K 



Ic 



a ,. a 

y correction j 



we may conclude that the maximum crack length has a 



interactive relation with the stress intensity factor K 



Itotal 



Following is the table which shows temperature mechanical properties for some materials used in 
engineering design . 



Table 1. 



Material 


a 


a 

ts 


^7. 


E 


AerMet 100 


1724 


1965 


126 


207 


Ti-6Al-AV 1 


869 


958 


87 


117 


Ti-6Al-AV2 


1007 


1034 


40 


130 


AISIA2,AQ 1 


1089 


1097 


110 


207 


AISIA3A0 2 


1476 


1896 


81 


207 


AISIAlAl 


945 


1062 


120 


207 


INCONEL 


1172 


1404 


96 


- 


18 M (250) 


1290 


1345 


176 


- 


2014 -T651 


455 


- 


24 


- 


2024 -T3 


345 


- 


44 


- 


7075 -T6 


572 


641 


24 


72 


NU,Fe,,P,,B,Si, 


800 


800 


12 


- 


SIC 


460 


460 


3.7 


72 



3. Construct algorithms for applications 

Here the author will present two procedures that are used to calculate the maximum pressure ( Pmax ) 
causing the fracture propagation and estimate the limitation of time while applying the continous or 
discrete internal pressure functions for the round metal vessel under the given constraints . 

3.1 Calculating the maximum pressure : 

The matrix A below is used to define the mechanical properties of each type of metal which is required 
for engineering design . 

pressure^^^ :=proc(a7, Bl, cl, dl, PI, typel, k) 



the data of problem will be entered into the procedure with the instructions below 

al : one-half crack depth (mm) , 

Bl : thickness of vessel cm) , 

cl : one-half crack length (mm) , 

dl : internal diameter(m) , 

PI : internal pressure (MPa) , 

typel : type of material ( inserted in the form of an element of the matrix A i.e A [k][i] ) 

k : the kth row ( see table 1 . ) . 



> restart; A ■=Matrix(\\ material, a ,a ,Kj^ ,E \ [AerMet J 00 , 1724 ,1965 , 

126 ,207 I [Ti6Al4VJ , S69 ,958 ,87 ,117 I [Ti6Al4V2 , 1007 ,1034 ,40 , 
130 ], [AISI43401 , 1089 , 1097 , 110 ,207 ], [AISI43402 , 1476 , 1896 , 81 , 
207 I [AISI4 147, 945 ,1062 ,120,207 ], [INC ONEL , 1172 ,1404 ,96,^^4 ], 
[18 M250, 1290 ,1345 ,176,7V:4 ], [2014 7657 , 455 , NA ,24 ,NA ], 
[2024 T3 , 345 ,NA ,44 ,NA ], [7075 T6 , 572 , 641 ,24 ,72 ], 
[Ni49Fe29P14B6Si2, SOO ,800 , \2 , NA ], [^/C,460 ,460 ,3.7,72 ]]) 

14 X 5 Matrix 

Data Type: anything 

Storage: rectangular 

Order: Fortran order 

> pressure^^^ :=proc(a7, Bl, cl, dl, PI, typel, k) local i, s, anew, cnew, Bnew, sigma , 

sigmal , alphal , alphax, P, Q, Ql, K , KIl , KIC ; with(pIots, implicitplot) : 

print ( " Start the program " ) ; print " The parameters are given : ",a = al,B = Bl,2c = 2 



(1) 



• cl, d = dl, P =P1, type = typel. 



a 



2-c 



evalf 



al 
2-cl ' 



3 , row = k 



nl /? / 

print { " Convert the units into meter : " ) ; anew ■= ; Bnew ■= ; cnew 



1000 



100 ' 



cl 



1000 



,print{a =anew, B =Bnew, c = cnew, d = dl);Kj-= alpha-P-y n-a 



printi "(1) The total stress intensity factor KI" = a[pha-P ^,^^^^^^^-J K-a ]; ; alphax 



1.12 



+ 1 



\ '^shape 



diameter 

2-B.., 

mickness 



+ 1 



; print 



' ^shape 



71 



4 I 2-c 



2\ 



pressure 

'(2) Correction factor a 
; sigma 



1.12 



Va 



d 



shape 



\ 2-B 



7 
33 



a 

\ y^ J 



pressure diameter 



2-B 



thickness 



print 



(3) The hoop stress a 



P -d 

pressure diameter 

9 • R 

thickness 



/»rz>7?r' (4) The shape factor Qshape " = Q^,\, print/ {"XnYos", 

" Read the data from the temperature mechanical properties table" ); print/ { "\n%s", " " ) ;; 
print(TYPE);fori from \ to 5 do print(s[i]=A[k][i]) ; s[i] ■= A[k][i] ; od; ; 
printi Material = s[ 1 ], Z ^ = 5[2 ], I^^ = s[3],KIc=s[4],E = s[5]y, ; ■,print/{ "\n%s", 

" Calculating the hoop stress sigma , the shape factor Q , correction factor alpha and the 
total stress intensity factor KI "); print/ {"\n%s" 



^11 II II 



); sigmal ■= evalf \ ^ , 5 



print {" From (3) "}; print 



(3) The hoop stress a 



pressure diameter 

2^. 



2-Bnew ' 



thickness 

print{" a "= sigma 1);; print/ {"XnVos"," ");printr 'From (4)"); 

2 . 

print [ " (4) The shape factor Q " = Q^^^^^); C5^^ -.= s[2]; Ql ■■= ^ -\j 



+ 



anew 
2-cnew 



2\ 



7 
33 



a 

\ ys J 



■ Ql := evalfiQl,!,)- -printi'' Qshape " =Q1)- ; 



f 



print/ i"\n%s'\ 



d 



diameter 



2-B 



+ 1 



thickness 



print [" From (2) ");print 
1.12 



'(2) Correction factor a 



1.12 



7a 



shape 



; alpha] 



dl 



, 2-Bnew 
"); print r From (\) 



+ 1 ; alpha 1 ■■= evalf {alpha 1, 3 ); 



print {" a " = alphal); ; print/ {"\n%s", 
printi "(1) The total stress intensity KI" = alpha.-P ^^^^^^^-^ n-a ]; X/7 ■= alphal-Pl 



pressure 

•7 Tt-awew ; print{KI = KI1);KI1 ■■= eval/{KIl, 4) ; print [KI = KI1); ; print/ ("\n%s'\ 

" ");/'r/>7//'("\n%s", ''!/• Compare the total stress intensity factor KI andKIc "); ; 

print/ {"\n%s'\ " ");X/C := s[4 ]; print {KIc=KIC);if KIl < KIC 

then print{ " The fracture will not occur / the crack propagation is stable , because KI < 

KIc " ); ; print/ rXnVos"," " );; print/ rXnVos"," "); print/ {"XnVos", 

" 2/. Determine the maximum pressure Pmax that causes the crack propagation " ) ; ; 

print/ {"\n%s"," ");print(" From (\) "); printi "(1) The critical stress intensity KIc" 

evalf ( Pf^^yix ' ^ ) ' Pf^f^f ( Pmax 



a 



[^]-^M..-V^^);^. 



MAX 



al-y] Ti- 



P 



MAX 



anew 



=Pmax ); ;;^n>7//-("\n%s", " ");/7r/>7//'("\n%s", 
" 3/. Consider the relation between Q and a " ) ; ; print/ { "\n%s", " " ) ; 
print(" From (2)"); print{ "Substitute x for alpha and y for Qshape into (2)"); 



print ( alpha = alphax ) ; x = evalf 



1.12 

V7 



dJ 



, 2-Bnew 



+ 1 ,4 



, print (%, plot 1) ,y = x-Pl 



■J n-anew;; printi %, plotl ) ;; print/ { "\n%s", " " ) ;, print/ { "\n%s", " " ) ; 
print{" Graphical relations between correction and shape factor (plotl) _ correction and 
stress intensity factor KI (plot2) " );;print{" End of program ") ; ; print/ ("\n%s"," "); 
^ 1.12 f dl 



plot 



eval/ 



. /x' 



, 2-Bnew 



+ 1 



:-PPj~K- 



anew 



, x = ..4,y = ..20, axes 



= boxed, labels = [ Q , Kf, alpha \ color = [blue, red], legend = [plotl, plot2], thickness 

= [2,1] 

else print( " Warning ! The fracture will occur / the crack propagation is unstable , 
because KI >= KIc " );,print{" End of program " );; fi; end: 



Example 3.1.1 

Consider a round vessel with technical data below : 

one-half crack depth a = 5.5mm ; thickness of vessel B = 20cm ; one-half crack length c = 4.5mm ; 

internal diameter d = 0.5m ; internal pressure P =362MPa ; type: AISI 4340 2 i.e A[6][l] 

Calculate : 

a. The stress intensity factor KI . 



b.The maximum static pressure that would cause the facture propagation . 

Compare the total stress intensity factor KI and KIc ( see table 1 . ) and give your conclusion 

Solution : 

Run the procedure * 



> pressure^JS.5, 20, 4.5, 0.5, 362,^[6][1 ], 6); 

" Start the program " 
" The parameters are given : ", a = 5.5, 5 = 20, 2 c = 9.0, t/ = 0.5, P = 362, type =AISI43402, 
1 a 



0.6\0,roM' = 6 

2 c 



Convert the units into meter 



a = 0.005500000000, B=^,c = 0.004500000000, d = 0.5 



" (1) The total stress intensity factor KI " = a/',,,.^^^,^^ J na 



pressure 



1.12 
'(2) Correction factor a " = 



^ 1 d 



U B+' 



/a 



shape 



•(3)Thehoops.ressa-i ^-— ''"— 



2 B.. , 

thickness 



1 2 

(4) The shape factor Qshape " = ~r ti 



3_ ]_ a^ 
442 



2 n2 2 



33 2 
a 

ys 



Read the data from the temperature mechanical properties table 

TYPE 



^1 


= AISI43402 




'2 = 


1476 




■^3 = 


1896 




^4 


= 81 




■^5 = 


= 207 



Material = AISI43402, Z = 1 476, Z = 1 896, X/c = 8 1 , £ = 207 

' ys ' ts ' ' 

Calculating the hoop stress sigma , the shape factor Q , 
correction factor alpha and the total stress intensity factor 
KI 

" From (3) " 



" (3) The hoop stress a 



1 pressure diameter 

~2 'b. 



thickness 



II _. II 



a "=452.50 



" From (4) 
" (4) The shape factor Q " = ^ 7t^ 



3 1a 

4 4 2 



7 a 



33 2 
a 

>'5 



" Qshape " = 3.10 



" From (2) " 
73.92 



1 diameter 



"(2) Correction factor a 



2 B 



+ 1 



thickness 



1089 71 
" a " = 1.43 



o 1 2 ^ 
3 1a' 

A A 2 



924 a 



VS 



" From (1) " 
(1) The total stress intensity ^^" = OiPp,.essurey '^^ 
KI = 3839069309 J~K 



pressure 



KI = 68.07 



1/ . Compare the total stress intensity factor KI and KIc 

KIc = Sl 
" The fracture will not occur / the crack propagation is stable , because KI < KIc 



2/. Determine the maximum pressure Pmax that causes the crack 
propagation 

" From (1) " 

" (1) The critical stress intensity KIc " = a Pj^f^^^J Tra 

Pwax = 430.90 



3/. Consider the relation between Q and &alpha; 



" From (2) " 
Substitute x for alpha and y for Qshape into (2) 



1.12 



1 diameter ^ 



2 5. 



a = 



thickness 



I 



2.520 , ^, 

vy 

y = 26.84663852 x J~K,pIot2 



" Graphical relations between correction and shape factor (plotl) _ correction and stress 
intensity factor KI (plot2) " 

" End of program " 



20-f 



a 10- 




O.KI 



plot! 



plot2 



Example 3.1.2 

Consider a round vessel with technical data below : 

one-half crack depth a = 4.5mm ; thickness of vessel B = 20cm ; one-half crack length c = 6mm 

internal diameter d = 0.8m ; internal pressure P = 560 MPa ; type: AISI 4340 1 i.e A[5][l] 

Calculate : 

a. The stress intensity factor KI . 

b.The maximum static pressure that would cause the facture propagation . 

Compare the total stress intensity factor KI and KIc ( see table 1 . ) and give your conclusion . 

Solution : 

Run the procedure pressure^^^{ , , , ) ; 



> pressure^J4.5, 20, 6, 0.8, 560,^[5][1 ], 5); 

" Start the program " 
" The parameters are given : ", a = 4.5, B = 20,2c= 12, t/= 0.8, P = 560, type =AISI4 3401, 
1 a 



2 c 



0.375, row = 5 



" Convert the units into meter : " 

a = 0.004500000000, B=—,c= ^— , d = O.S 

5 ' 500 ' 



(1) The total stress intensity factor KI " = OiP^^,^^^^^^^ J na 

1.12 



^ 1 d 



"(2) Correction factor a 

" (3) The hoop stress a 
" (4) The shape factor Qshape " = 



2 B 



pressure 
+ 1 



y '--shape 
1 pressure diameter 



B. 



thickness 



1 2 



4 4 ^2 



2 N 



7 a 



33 2 
a 

V5 



Read the data from the temperature mechanical properties table 

TYPE 



h 


= AISI43401 




'2 = 


= 1089 




'h = 


= 1097 




S4- 


= 110 




'5- 


= 207 



Material =AISI4 3401, E = 1089, Z = 1097, X/c = 1 10, £ = 207 

' ys ^ ts ' ' 



Calculating the hoop stress sigma , the shape factor Q , 
correction factor alpha and the total stress intensity factor 
KI 



" From (3) 



(3) The hoop stress a 



1 pressure diameter 

2" B.., 

thicliness 



" From (4) 

1 2 
(4) The shape factor Q " = — k 



442 



2 N 



7 o_ 

33 2 



VS 



" Qshape " = 1.73 



" From (2) " 
73.92 



"(2) Correction factor a 



1 diameter ^ 



thickness 



1089 ^ 



^1,1 (^ f 924 (5 



2 

a 

ys 



a " = 2.55 



" From (1) " 
" (1) The total stress intensity KI " = aP^^.g^^^^^g V Tia 

X/ = 95.79315215 7^ 



pressure 



Ki=\m.% 



1/ . Compare the total stress intensity factor KI and KIc 

nc = 110 

" Warning ! The fracture will occur / the crack propagation is unstable , because KI >= KIc 

" End of program " (2) 



3.2 Estimating the interval of time while applying the internal dynamic pressure 

For the next procedure pressure^^^^ :=proc(a7, Bl, cl, dl, PI, typel, k, dim, tO, yO, yl ) 

we enter the data of the problem with following instructions 

al : one-half crack depth (mm) , 

Bl : thickness of vessel cm) , 

cl : one-half crack length (mm) , 

dl : internal diameter(m) , 

PI : internal pressure function ( continous /discrete ) (MPa) , 

typel : type of material ( inserted in the form of an element of the matrix A i.e A [k][i] ) , 

k : the kth row ( see table 1 . ) , 

tlim : the upper bound of interval ( , tlim ) that is a constraint of time used to solve the inequality KI 

< KIc, 

to, yO , yl : the minimum and maximum values of t and y used to plot the graph of funtions KI and 

the polynomial KITaylor . 



> pressure^^^^ :=proc (a7, Bl, cl, dl, PI, typel, k, tlim, tO, yO, yl ) global P2 , KIC, KIl, KI2 
local St, st2, ineq, ineq2, i, s, anew, cnew, Bnew, sigma , sigmal , alphal , alphax, Q, P, 



Q1,K ,with[plots,implicitplot) :prmt{" Start the program " ); 

print \ " Parameters of the metal vessel and internal pressure P : ", a=aJ, B = BJ, 2 c = 2 



ci ( cil 

■ cl,d = dl,P =Pl,type = typel, — — =evalf\ — — r, 3 \,row = k 



2-c 



2-cl ' 



n t J^ I 

print { " Convert the units into meter : " ) ; anew ■■= ; Bnew ■= ; cnew 



1000 



100 ' 



cl 



,printia=anew,B = Bnew, c = cnew, d = dI),P2 ■= Pr,Kj-= alpha-P ■ ^ n-a ; 
printi "(1) The total stress intensity factor KI " =alpha-P -^ n-a ); ^alphax 



1.12 



Va 



shape 



diameter 

2-B.., 

thickness 



pressure 



d 



diameter 



2-B 



+ 1 



thickness 



■.Q--= 



+ 1 



71 



print 



(3) The hoop stress a 



; print 


"(2) Correction factor a 


\ 

3 f a \ \ 


2 

7 


' a ^ 


U I 2-c J J 33 


a 


pressure diameter 




9 • /? 

thickness 


/ 


5 5 





1.12 



Va 



shape 



Sigma := 



2-5' 



/7rz>7^( " (4) The shape factor Q " = Q^^^p^); print/ ( "\n%s", 

" Read the data from the temperature mechanical properties table" ) ; print/ { "\n%s", " " ) ;; 
print(TYPE);fori from \ to 5 do print{s[i]=A[k][i]) ; s[i] ■■= A[k][i] ; od; ;; 
printi Material = s[ 1 ], Z ^ = ^[2 ], Z^^ = s[3], KIc = s[4], E = s[5]y, ;,print/{ "\n%s", 

" Calculating the hoop stress sigma , the shape factor Q , correction factor alpha and the 
total stress intensity factor KI "); print/ {"\n%s" 



^11 II II 



); sigmal ■= eval/\ ^ , 5 



print {" From (3) "}; print 



(3) The hoop stress a 



pressure diameter 

2^. 



2-Bnew ' 



thickness 

printi" cr "= sigmal);; print/ {"XnVos"," "); print (" From (4)"); 

2 . 

print [" (4) The shape factor Q "=Q,,,p,); o,^ == ^[2]; QJ ■■= ^ ■ 

ol 



4 



+ 



anew 



2\ 



7 
33 



/ _, x2 



a 



f 



2 -cnew 
print/ { "\n%s", " " ) ; print ( " From (2) " ) ; print 

diameter , ^ j » , 1.12 ( dJ 



;Ql:=eval/{Ql,3);;printr Q " =Q1);; 



'(2) Correction factor a 



1.12 



d 



Va 



shape 



2-B 



+ 1 



thickness 



; alpha] 



Tqt 



2 -Bnew 



+ 1 ; alpha 1 ■= evalf {alpha 1, 3 ); 



print ( " a " = alphal ) ; ; print/ { "\n%s", " " ) ; print ( " From ( 1 ) " ) ; 

print\ " (1) The total stress intensity factor KI " =alpha-P-7 Ti-a ); KIl ■= alpha] -P J 

•7 n-anew ; print(KI = KI1);KI1 ■■= eval/{KIl, 4) ; print {KI = KII); ; print/ {"\n%s", 
" ");/'r/>7//'("\n%s", " 1/. Compare the total stress intensity factor KI and KIc "); ; 

print/ r\n%s", " ");KIC ■= s[4];print{KIc=KIC); 

print( " The fracture will not occur / the crack propagation is stable , as KI < KIc " ) ; 
ineq ■■= KIl < KIC ;print{" Analytical inequality ", ineq);KI2 ■= convert[taylor{KIl, t 



= ^0),/?o/v«om);/?n>7^(" Taylor polynomial expression of KI at t = " , tO, KITaylor 

= KI2); ineq2 '■— convert(KI2, polynom) < KIC; 

print { " The analytical inequality obtained by using Taylor polynomial expression ", 
meq2); st2 ■■= solve { {meq2, t > 0, t < dim}, t) ; 

print {" Estimating time-limit based on Taylor polynomial expression of KI ", st2); 
print/ ( "\n%s", " ");st ■■= solve { {ineq, t > 0, t < tlim}, t) ; 

/>n>7?(" Estimating time— limit based on analytical method for KI ", st); print/ {"\n%s", 
" "); print/ r\n%s", 
" 2/. Determine the maximum pressure Pmax that causes the crack propagation " ) ; ; 

print/ {"\n%s'\" ");print(" From (\) ");print( "(1) The critical stress intensity KIc " 



= o^[^]-^M«x-V^^);^A^ == . '^Pa4AX '■= ^^^V{Pmax^ 5); print[Pmax 

al -^ Ti-anew 

=Pmax ); ;/^n«?/'("\n%s", " ");/^n>7r/("\n%s", 
" 3/. Consider the relation between Q and a " ) ; ; print/ [ "\n%s", " " ) ; 
print(" 'From (2)"); print { "Substitute x for alpha and y for Q into (2)"); 



print ( alpha = alphax ) ; x = evalf 



1.12 



dl 



+ 1 



; print { %) ;; print/[ "\n%s" 



, 2-Bnew 

" ");■ print/ {"\nyos'\" ");''"' 

print( " Graphical relations between t and stress intensity factor (plotl) _ Taylor series 
expansion of stress intensity factor KI (plot2) " ).,;print{" End of the program " );;;; 
print/{ "\n%s", " " );plot( [KIJ , KI2 lt = ..tlim, y =y0.yl, axes = boxed, labels =Vt\ 
KI '], color = {blue, green], legend= [Klplotl, KITaylorplot2], thickness= [3, 3 ], linestyle 
= [solid, dot]); end: 

> P^^ssure^j^^j^^ ■= proc ( tKcor, ymin, ymax ) local pK ; 

print( " The value of KIc is given from the engineering design requirement ( see table 1. ) : 
KIc =",KIC); 

print{ " The expresion of continous or discrete form of KI and KITaylor " ) ; 
print[ " KI = ", KI[t)=KIl); print( " KITaylor = ", KITaylor {t)=KI2); 

; plot{ [KII, KI2, KIC ], t = .. tKcor, y =ymin . .ymax, color = [ blue, 
green, red], thickness = [3, 3, 3], legend = [ ^bhieKIV, greenKITaylor\ ^edKIc^], linestyle 
= [solid, dot, dash]); end: 

> pressiirCpi^^.^ ^^ •■=proc{tcor, ymin, ymax) local pp ; 

print(" The expresion of continous or discrete form of internal pressure " ); 
print[" Pressure function = ",/(0 ^P2); ;plot{P2, t = .. tcor,y 

=ymin . .ymax, color = brown, thickness = 3, legend= ^brownP ') ; end : 



Example 3.2.1 

Consider a round vessel with technical data below : 

one-half crack depth a = 5mm ; thickness of vessel B = 25cm ; one-half crack length c = 8mm ; 

internal diameter d = 0.8m ; internal pressure P = 

16.13084112 ln(r + 15.22) +2.5 e^'"(°°^^°^4^i°3i ' + °-^^428060i2) -9600.268 ^^^ . ^^^ . Ti6A14V2 i.e 

A[4][l] 

a. Find the expression of stress intensity factor KI and Taylor polynomial of KI . 

b. Evaluating the time-limit for applying the continous / discrete form of internal dynamic pressure . 



Solution : 

Run the procedure P^essure^^^^ (,,,); 

> pressurej,^^^{5, 25,8,0.8, 16.13084112 ln(r + 15.22) 

^2 5 esm(0.06507451031/ + 0.6642806012) - 9600.268^ ^ ^^ ^ ^ ^ ^^ ^^ ^qqqq^ 5qqq^ q^ 3q ) . 

" Start the program " 
" Parameters of the metal vessel and internal pressure P : ", a = 5, B = 25, 2 c = 16, d = O.S, P 



sin(0.06507451031 / + 0.6642806012) - 9600.268 



16.13084112 ln(r + 15.22) +2.5 e' 

Ti6Al4V2, — — = 0312, row = A 

2 c 

" Convert the units into meter : " 

a=^^,B=—,c=-^,d=0.% 
200 ' 4 ' 125 ' 

" (1) The total stress intensity factor KI " = oc/" ^^.^.^^^.^ ^ na 



type 



pressure 



1.12 



1 diameter 



"(2) Correction factor a 



" (3) The hoop stress a 



2 B 



+ 1 



thickness 



iQs^ 



shape 
1 Ppressuredd^ameter 



B 



(4) The shape factor Q 



1 2 
7^ 



f 3_ J_ a^ 
4 4 2 



thickness 
2 \2 



7 a 



2 \ 



33 2 
a 

_V5 y 



shape 



Read the data from the temperature mechanical properties table 

TYPE 



•^1 


= Ti6Al4V2 




^2 = 1007 




^3 = 1034 




^4 = 40 




5g=130 



Material= Ti6Al4V2, Z = 1007, 2 = 1034, KIc = 40, £ = 130 

' ys ^ ts ' ' 



Calculating the hoop stress sigma , the shape factor Q , 
correction factor alpha and the total stress intensity factor 
KI 



" From (3) 



" (3) The hoop stress a 



1 pressure diameter 



fi ,— II 



a " = 25.8101n(r + 15.22) +4.0000e' 



thickness 
sin(0.065075 / + 0.66428) - 9600.3 



" (4) The shape factor Q 



" From (4) 



1 2 



2 \ 



4 4 


a 


7 


o 


c'] 


33 


2 

a 








j^i ; 



shape 



Q " = 1.77 -2.09 10"' (25.8 ln(? + 15.2) +4 qo e""(°°'''^ + °-''"' -''°°)' 



'(2) Correction factor a 



a 



" From (2) " 
1.12 



1 diameter ^ 



2 5. 



thickness 



J 




1.77 -2.09 10"^ (25.8 ln(r + 15.2) +4.00e' 



sin(0.0651/ + 0.664) -9600.^ 



" From (1) 



" (1) The total stress intensity factor KI" = aP ^ na 
KI= (o.l455000000 ( 16.13084112 ln(? + 15.22) 

,r^r sin(0.06507451031 / + 0.6642806012) -9600.268^ ry j W 

1.77 -2.09 10-' (25.8 ln(r + 15.2) +4.00 e^'"(0°6^i^ + 0-664) -9600.)2 

0.3648 (16.13 ln(r+ 15.22) + 2.5 e""(°°6507^ + 06643) -9600.) 



KI 



1.77 -2.09 10"' (25.8 ln(r + 15.2) +4.00e' 



sin(0.0651/ + 0.664) -9600. 



1/. Compare the total stress intensity factor KI and KIc 

KIc = 40 
" The fracture will not occur / the crack propagation is stable , as KI < KIc " 
Analytical inequality ", 



0.3648 ( 16.13 ln(r + 15.22) +2.5 e""'°°^^°^^ + °'^'^^^) -'^^°°) ^ ^^ 

7 1.77 -2.09 10-' (25.8 ln(r+ 15.2) +4.00 e^'"(°°''i^ + °-^^4'-^^°°f 
" Taylor polynomial expression of KI at t = ", 5000, KITaylor = 33.34426980 

+ 0.0008894882453^-8.849970677 10"^ (^ - 5000)^ + 1.175367970 10^^^ (^-5000)^ 

-1.756728539 10^^^ (^ - 5000)^ +2.801127987 10^^^ (? - 5000)^ 
" The analytical inequality obtained by using Taylor polynomial expression ", 

0.0008894882453 t - 8.849970677 10^^ (t - 5000)^ + 1.175367970 10^" {t - SOOOf 

-1.756728539 10^^^ (^ - 5000)"* +2.801127987 10^^^(^-5000)^ < 6.65573020 
" Estimating time-limit based on Taylor polynomial expression of KI ", {0. < t, t 
< 8155.424162} 



Estimating time-limit based on analytical method for KI 



2/. Determine the maximum pressure Pmax that causes the crack 
propagation 

" From (1) " 
" (1) The critical stress intensity KIc " = a Pj^f^_^^J t^ci 



Pmax =109. 61 J \. 11 -2.09 10"' (25.8 ln(r + 15.2) +4 qo e""(°°'''' + °-'''^ -''°°)' 



3/. Consider the relation between Q and &alpha; 

" From (2) " 
" Substitute x for alpha and y for Q into (2) 



1.12 
a= 



1 diameter .. 

2" B.., 

thickness 



-y '-■shape 

^ 2.912 

' V7 



" Graphical relations between t and stress intensity factor (plotl) _ Taylor series expansion of 
stress intensity factor KI (plot2) " 

" End of the program " 



80H 



70- 



60- 



50- 



KI 40- 




10000 



Klplotl ' ' ' KITaylorplot2 



" The value of KIc is given from the engineering design requirement ( see table 1. ) : KIc 
40 

" The expresion of continous or discrete form of KI and KITaylor " 

0.3648 (16.13 ln(r + 15.22) +2.5 e^'"(°°^^°'' + °-^^4^) -^^°°) 



"KI=",Z7(0 



1.77-2.09 10"^ (25.8 ln(? + 15.2) +4.00e' 



sin(0.0651/ + 0.664) -9600. 



r 



KlT?iy\or=\ KITaylor [t) =33.34426980 +0.0008894882453 ^ - 8.849970677 lO"'' [t 
- 5000)^ + 1.175367970 10"" [t - 5000)^ - 1.756728539 10"^^ {t - 5000)^ 
+ 2.801127987 10"^^ (^ - 5000)^ 



80 n 



60- 



y 40 








20-' 



-1 1 1 r 



-1 1 1 r 



5000 



10000 



-1 1 1 r 



15000 



blueKIl 



greenKITaylor 



redKIc 



> pressurep,^^.„^J\000, 00, 122.28) ; 

" The expresion of continous or discrete form of internal pressure 
" Pressure function = \f{t) = 16.13084112 \n(t + 15.22) 

,nr sin(0.06507451031 1 + 0.6642806012) - 9600.268 



120- 



100- 



80- 



y 60- 



40- 



20- 



0' 







200 



400 



600 



800 



1000 



t 



brownP 



Example 3.2.2 

Consider a round vessel with technical data below : 

one-half crack depth a = 5mm ; thickness of vessel B = 25cm ; one-half crack length c = 8mm ; 
internal diameter d = 0.8m ; internal pressure P is given from the data table below ; type: 18Ni 
(250) i.e A[9][l] 



18Ni 
(250) 


adjust 
coeff 


1/ 

948\ 
.89\ 

2 




















t 





5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


MPa 





150 


170 


140 


165 


220 


340 


800 


920 


900 


820 



a. Find the expression of stress intensity factor KI and Taylor polynomial of KI . 

b. Evaluating the time-limit for applying the discrete form of internal dynamic pressure P 



Solution 



Fitting the data points by 

jhumPower ■■= PolynomialInterpolation{[[0, 0.0], [5, 150], [10, 170], [15, 140], [20, 165], [25, 
220], [30,340], [35,800], [40,920], [45 ,900 ], [50 ,820 ]], tJorm= power); 



> with[CurveFittmg) :with(pIots) .JhumPower -■= Polynomiallnterpolation([[0, 0.0], [5, 

150], [10, 170], [15, 140], [20, 165], [25, 220], [30, 340], [35, 800], [40, 920], [45 
900 ], [50 ,820 ]], t, form = power); 

JhumPower := -6.439506133 10"^°^^° + 1.544550255 10"^ ^^ - 0.00001579195757 /^ 
+ 0.0008993206290 t^ - 0.03 12803 1829 f + 0.6852472173 t^ - 9.393835424 t^ 
+ 77.20204438 P - 344.1569100 f^ + 651.5571380 t 

> pointploti {[0,0.0], [5, 150], [10, 170], [15, 140], [20, 165], [25,220], [30,340], [35, 

800], [40,920], [45 ,900 ], [50 ,820 ]}, axes = framed); 

900-1 * • 



800- 



700- 



600- 



500- 



400- 



300- 



200- 



100- 



0-lf^ 



10 



20 



30 



40 



50 



(3) 



Kunthis procedure pressiire^^^^^{, , , ) with the approximate polynomial friumPower 



> pressure^^^^ 



5, 25, 8, 0.8, -^'^^l^^' ,A[9][\l 9, 80, 47,-50, 100 



5 ■— -- 5 "--J " • "--J 



" Start the program " 
" Parameters of the metal vessel and internal pressure P : ", a = 5, B = 25, 2 c = 16, d = O.S, P 

-6.786342525 10"^^?^° + 1.627740833 10"^°?^ - 1.664252366 10" V 

+ 9.477586793 10"^ t^ - 0.00003296509854 f + 0.0007221551210 t^ 

- 0.009899794098 t^ + 0.08136020154 ^ - 0.3626934466 f^ + 0.6866504700 t, type 

= \SNi250, — — = 0.3 12, row = 9 

2 c 

" Convert the units into meter : " 

a=^— ,5=— ,c=^— ,£/=0.8 
200 ' 4 ' 125 ' 

" (1) The total stress intensity factor KI " = ocP ^^^^^^^ ^ na 



1.12 



"(2) Correction factor a 



(3) The hoop stress a 



( 1 d 

1 diameter 



2 B 



+ 1 



thickness 



/e 



■shape 



1 pressure diameter 



(4) The shape factor Q 



1 2 
7^ 



4 4 ^2 



thickness 
2 \ 



7 a 



2 \ 



33 2 
o 

y ' shape 



Read the data from the temperature mechanical properties table 

TYPE 



s^ = \%Ni250 


'2 = 


= 1290 


^3 = 


= 1345 


^4 


= 176 


•^5 


=NA 



Material= 18 Ni250, E = 1290, Z = 1345, KIc = \16,E = NA 

' ys ' ts ' ' 



Calculating the hoop stress sigma , the shape factor Q , 
correction factor alpha and the total stress intensity factor 
KI 



" From (3) 



(3) The hoop stress a 



1 pressure diameter 



thickness 



II _. II . 



1.0858 10"^^?^° +2.6043 10"^°?^ -2.6629 10"^ r^ + 0.0000015164 /"^ 



0.000052745 i' + 0.001 1554 f' - 0.015840 t + 0. 13018 ? - 0.58030 ? + 1.0986 t 



" From (4) 



" (4) The shape factor Q 



1 2 
7^ 



2 N 



4 4 2 



7 o 



2 N 



33 2 
a 

>'5 ; 



shape 



Q " = 1.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^° ^^ -2.66 10"V + 0.00000152 ^^ 
- 0.0000527 i" + 0.00116 i - 0.0158 ?^ + 0.130 ? - 0.580 ^^ + 1.10 



" From (2) " 
1.12 



1 diameter 



"(2) Correction factor a 



2 B.. , 

thickness 



+ 1 



1 2 



4 4 ,2 ^ 



7 a 



2 \ 



33 2 
a 

ys ) 



shape 



" a " = 2.91/ 

(l.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^°?^ -2.66 10"^?^ 

+ 0.00000152 ? - 0.0000527 i' + 0.00116 ?^ - 0.0158 i" + 0.130 ? - 0.580 /^ + 1.10 ?) ) 

1/2 



" From (1) 



" (1) The total stress intensity factor KI " = aP ^ na 
KI= (o.l455000000 (-6.786342525 10"^^^^° + 1.627740833 10"^° ^^ - 1.664252366 \0'^ t^ 
+ 9.477586793 10"^ P - 0.00003296509854 f + 0.0007221551210 t^ 
- 0.009899794098 t^ + 0.08136020154 P - 0.3626934466 f^ + 0.6866504700 t) 

(l.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^°?^ -2.66 \0'^ t^ 

+ 0.00000152 f - 0.0000527 f + 0.001 16 /^ - 0.0158 t^ + 0.130 P - 0.580 P + l.lOt) ) 

ill 

Z7= (0.3648 (-6.786 10"^"^^^° + 1.628 10"^° ^^ - 1.664 10"V + 9.478 10"^ ?^ - 0.00003297 Z' 
+ 0.0007222 r' - 0.009900 / + 0.08136 P - 0.3627 ? + 0.6867 t) ) / 



(l.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^°?^ -2.66 \0'^ t^ 

+ 0.00000152 f - 0.0000527 f + 0.001 16 P - 0.0158 t^ + 0.130 P - 0.580 ;^ + 1.10 r) ) 

1/2 



1/. Compare the total stress intensity factor KI and KIc 

KIc = 176 
" The fracture will not occur / the crack propagation is stable , as KI < KIc " 
" Analytical inequality ", (0.3648 ( -6.786 10"^^^^° + 1.628 10"^° ^^ - 1.664 10"^ ^^ 

+ 9.478 10"^ t^ - 0.00003297 f + 0.0007222 t^ - 0.009900 t^ + 0.08136 f^ - 0.3627 ^ 

+ 0.6867 t)) / 

(l.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^°?^ -2.66 \0'^ t^ 

+ 0.00000152 P - 0.0000527 f + 0.00116 f' - 0.0158 t^ + 0.130 P - 0.580 P + l.lOt) ) 

ill 

< 176 

" Taylor polynomial expression of KI at t = ", 47, KITaylor = -2377.576003 

+ 51.69382641 ^ + 63.00700353 (? -47)^ + 88.08829310 (^-47)^ 

+ 130.3049045 (? -47)^ + 198.6925598 (^-47)^ 
" The analytical inequality obtained by using Taylor polynomial expression ", 51.69382641 t 

+ 63.00700353 (^ -47)^ + 88.08829310 (r-47)V 130.3049045 (^-47)'* 

+ 198.6925598 (^-47)^ < 2553.576003 
" Estimating time-limit based on Taylor polynomial expression of KI ", {0. <t,t 

< 47.63775925} 



" Estimating time-limit based on analytical method for KI ", {0. < r, ?< 47.52889066}, 
{71.07919106 <t,t < 72.30034108} 



2/. Determine the maximum pressure Pmax that causes the crack 
propagation 

" From (1) " 
" (1) The critical stress intensity KIc " = a P^,f^^^J T^a 
Pmax 

= 482.55 (l.77 -1.27 10"^ (-1.09 10"^^/° +2.60 10"^° r^ -2.66 10"V 



9 >v 

+ 0.00000152 t^ - 0.0000527 i' + 0.00116 i - 0.0158 t^ + 0.130 f' - 0.580 ? + l.lOt) ) 

1/2 



3/. Consider the relation between Q and &alpha; 



" From (2) " 
" Substitute x for alpha and y for Q into (2) 



1.12 



a = 



1 diameter -. 

thickness j 



_ 2.912 

' V7 



Graphical relations between t and stress intensity factor (plotl) _ Taylor series expansion of 
stress intensity factor KI (plot2) " 

" End of the program " 



loo-r 



KI 




Klplotl • ■ KITaylorplot2 



> pressure ^^^jj^^ ( 80, - 50, 1 90 ) 

" The value of KIc is given from the engineering design requirement ( see table 1. ) : KIc = ", 
176 

" The expresion of continous or discrete form of KI and KITaylor " 

"KI=",X/(0 = (0.3648 (-6.786 10"^-^^^° + 1.628 10"^° ^^ - 1.664 10"^ r^ + 9.478 \Qi\^ 

- 0.00003297 f + 0.0007222 ^ - 0.009900 i" + 0.08136 f^ - 0.3627 ? + 0.6867 t) ) / 

(l.77 -1.27 10"^ (-1.09 10"^^^^° +2.60 10"^° ^^ -2.66 10" V 

+ 0.00000152 t^ - 0.0000527 /' + 0.00116 i - 0.0158 t^ + 0.130 r' - 0.580 ? + l.lOt) ) 

1/2 



" KITaylor =",X/ray/or(0 = -41.31589512 + 1.146168336^ + 0.1478373550 [t-AOy 
+ 0.01574740070 (? -40) V 0.001656396686 (? -40)^ + 0.0001763158242 (?-40)^ 



150- 



y 



100- 




blueKIl 



greenKITaylor 



redKIc 



> pressurepj^^^^^ JlO,-05, I) -pointploti {[0,0.01 [5, 150], [10, 170], [15, 140], [20, 165], 

[25,220], [30,340], [35,800], [40,920], [45 ,900 ], [50 ,820 ]}, axes = framed); 
" The expresion of continous or discrete form of internal pressure " 

" Pressure function = ",/(0 = -6.786342525 10"^^?^° + 1.627740833 10"^° ^^ 

- 1.664252366 10"^ ? + 9.477586793 10"^ t^ - 0.00003296509854 /^ 

+ 0.0007221551210 f' - 0.009899794098 f + 0.08136020154 f^ - 0.3626934466 ? 
+ 0.6866504700 r 



y 0.5- 




■0.5 -■ 



brownP 



900- 



800- 



700- 



600- 



500- 



400- 



300- 



200- 



100- 



0^^ 



10 



20 



30 
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Notice that if we use Taylor series expansion of the approximate function for internal dynamic pressure 
at t = 47 it gives the same results of the time4imit calculated by analytical and series methods . ( 
0. <t,t < 47.63775925 , 0. <tj < 47.52889066 ) 



Conclusion 

Through the examples presented above, we find that the Maple tools has a lot of utilities in engineering 
computation and estimating the technical specifications . Note that in the case of applying internal 
dynamics pressure of discrete data we can use some softwares to find the best fitted functions , then get 
the expressions of the approximate functions to enter the procedures . 

For the example 3.2.2 above, using the software Curve Expert , after entering the data we find the 
approximate expressions with parameters such as Standard Error, Correlation Coefficient ... From which 
we can choose the appropriate functions for the next calculation . 
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Fig 1. Data plot . 



1. Logistic Model: y=a/(l+b*exp(-cx)) 

Standard Error: 1 14.9209641 

Correlation Coefficient: 0.9580550 

Comments: 

The fit converged to a tolerance of le-006 in 36 iterations. No weighting used. 



S- 114.92096413 

r= 0.95805603 
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Fig 2 . Graph of the logistic fit . 



2. The 4th Degree Polynomial Fit: y=a+bx+cx^2+dx^3+exM 

Coefficient Data: 

a= -0.59440559 

b= 57.013598 

c= -6.1678322 

d= 0.23554002 

e = -0.0025734266 



Standard Error: 75.4632160 

Correlation Coefficient: 0.9866303 

Comments: 

Linear regression completed successfully. No weighting used. 



S=75l 46321600 
r =0.98663002 
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Fig 3 . Graph of the 4th Degree Polynomial Fit . 



3. Lagrangian Interpolation: y=a+bx+cx'^2+dx^3+exM+fx^5+gx^6+hx^7+ix^8+jx^9+kx^l0 
Coefficient Data: 



a = 
b = 
c = 
d = 
e = 
f= 







651.55714 

-344.15691 

77.202045 

-9.3938355 

0.68524722 

-0.031280319 

0.00089932063 

-1.5791958e-005 

1.5445503e-007 
k= -6.4395062e-010 

Standard Error: 0.0000000 
Correlation Coefficient: 1.0000000 
Comments: 

Interpolation completed successfully. Since interpolations pass through each data point, the error is 
and the corr. coeff is 1. 
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Fig 4 . Graph of the Lagrangian Interpolation fit . 
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